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Abstract. - A Ginzburg-Landau theory is presented on surfactants in polar binary mixtures, 
which aggregate at an interface due to the amphiphilic interaction. They can be ionic surfactants 
coexisting with counterions. Including the solvation and image interactions and accounting for a 
finite volume fraction of the surfactant, we obtain their distributions and the electric potential 
around an interface in equilibrium. The surface tension is also calculated. The distribution of the 
adsorbed ionic surfactant is narrower than that of the counterions. The adsorption is marked for 
hydrophilic and hydrophobic pairs of ionic surfactant and counterions. 



Introduction. — Surfactant molecules can be trapped 
at an interface, reducing the surface tension and giving 
rise to mesoscopic structures, as their density in the bulk 
is increased [1,2]. Oil- water interfaces containing surfac- 
tants are usually treated as infinitesimally thin surfaces 
characterized by the surface tension 7 dependent on the 
adsorbed surfactant density F. In this letter, we present 
a diffuse-interface model or a Ginzburg-Landau theory. 
Some authors presented Ginzburg-Landau free energies 
which include higher-order gradients of the composition or 
the surfactant density to describe dynamics of mesophase 
formation [3-6]. However, the gradient expansion in the 
free energy breaks down when the surfactants are strongly 
adsorbed onto an interface or when the trapping energy ea 
of a surfactant molecule much exceeds the thermal energy 
T. We will propose a model in which surfactant molecules 
are treated as rods and their two ends can stay in very dif- 
ferent environments. In this letter, we examine only the 
equilibrium properties following from our model. 

Moreover, surfactants can be ions or have charges in 
many important systems in biology and technology. The 
adsorption of surfactant and counterions in such cases has 
not yet been well studied. Hence we will treat ionic sur- 
factants and counterions in binary mixtures, grafting the 
amphiphilic interaction onto our previous theory of elec- 
trolytes which accounts for the solvation effect [7, 8] . 

We may define the excess adsorption F generally for 
doped particles such as surfactants and ions. For a diffuse- 
interface, we determine the interface position in terms of 
the order parameter ^l^{z) (the composition difference in 



this work) as 2:1,1 = /q dz[^(z) — '0^]/A?/' in a finite system 
in the region Q < z < L, where all the quantities change 
along the z axis [9]. The 4'{z) takes the two bulk values 
■0^ and ipi3 with A?/; = ipa — "4^(3 > on the two sides of 
the interface, while the bulk values of the doped-particle 
density n{z) are denoted by na and np. The phase a is 
more polar or water-rich, while the phase (3 is less polar 
or oil-rich. Then, 



T ^ I dz 
'0 



dz[n{z) — na\ + I dz[n(^z) — np], (1) 



where An = Ua — np and the integrands tend to zero far 
from the interface. At very small n, the surface tension 
decreases as 7 = 70 — TV, where 70 is the surface tension 
without doping [2,9]. For electrolyte systems [8] there is 
also an electrostatic contribution. 



7 ^ 70 - TF - / dze|V$| 
87r J 



(2) 



where e is the dielectric constant and $ is the electric 
potential tending to constants far from the interface. We 
will derive eq. (2) including the amphiphilic interaction. 
The Boltzmann constant will be set equal to unity. 

Nonionic surfactant in a binary mixture. We 

first consider a polar binary mixture (water and oil) with 
a small amount of neutral surfactant. The water, oil, 
and surfactant densities are ua, ns, and ni, respectively. 
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Assuming the same molecular size a for water and oil, 
their volume fractions arc 6 a = cl'^iT'A-, <t>B = o?nB, and 
A^'iCi = V\ni, respectively, where v\ is the surfactant 
molecular volume. We define the normalized surfactant 
density ci = a?ni. Here the volume ratio A^i — vi/a^ 
can be arbitrary. We assume the space-filling condition 
(pA + <l>B + Nici = 1; then, 

ct>A = {l-NiCi)/2 + i^, (/.B = (l-A^iCi)/2-V;. (3) 

Here 2ip = (j)A — (j^B the composition difference of the 
solvent. We may set (j)A — 1/2 + ip and <pB — 1/2 — 
ip neglecting the surfactant volume fraction for iViCi ^ 
1/2 — The free energy F is the space integral of its 
density /non of the form. 



J non 



4>A In (pA 
+Ci In ci 



'/>B In (/ 



J + X<t>A(t>B 

Ci InZam, 



(4) 



where x is the temperature-dependent coefficient, the co- 
efficient C is positive, and gi represents the interaction be- 
tween the surfactant and the composition difference. The 
last term represents the amphiphilic interaction, where 
Za.m is the partition function of a rod-like dipole with its 
center at the position r. We assume that the surfactant 
molecules take a rod-like shape with length 2(. consider- 
ably longer than a. It is given by the following integral on 
the surface of a sphere with radius i, 



,M = / 



4tt 



exp 



— Watpir + iuj+Watpir — iu) 



, (5) 



where u is the unit vector along the rod direction and 
/ dfl represents the integration over the angles of u. The 
two ends of the rod are at r + £u and r — £u under the 
influence of the solvation potentials given by Tw.^ijj{r + 
iu) and —T'Wi^'iJj{r — £u), respectively. It is instructive 
to examine the case in which ijj{r) varies slowly. That 
is, if the expansion ip{r + lu) — ■ip{r — £u) = 2£u ■ Vip is 
used, the last term in cq. (4) becomes — Iwi^f^jV-iApCi,. 
This gradient form was assumed in the literature [3,4]. 
Together with the fourth term in eq. (4), the coefficient 
in front of |VV'P vanishes for m = til (a Lifshitz point), 
where 

nL = SC/[4a''i£waf]. (6) 

For m > riL, a homogeneous solution is unstable at a 
finite wave number, leading to a mesophase. If C a^, 
we have a^riL ~ {a/£wa)^, so til is small for ^ a- 

The above gradient expansion cannot be used around 
an interface far from the critical point or for ^ < 2£. If 
a surfactant molecule is trapped at an interface with its 
hydrophilic (hydrophobic) end in the water-rich (oil-rich) 
region, the free energy decreases by = Tw^Aip. Strong 
adsorption occurs for large ea./T = Wa^ip. In the one- 
dimensional (ID) case, where all the quantities vary along 
the z axis, Z^m becomes 



dC 

2£ 



exp[-Watp{z + Q + Wa^jiz - Q], (7) 



where we have replaced tp{r ± £u) in eq. (5) by 'il){z ± Q) 
with C = £uz- In the thin interface limit ^ <C we place 
the interface at z = to find 

^am(^) = 1 + (1 - |^|/£)[cOsh(«;aA^) - 1] (8) 

for \z\ < while .^am = 1 for \z\ > £. Furthermore, in the 
dilute limit NiCi <C 1 — we have ci{z) = C\aZan\{z) 
for 2; < and c\{z) = ciisZamiz) for z > 0, where cia and 
ci/3 = e-fi^'^cia are the bulk surfactant densities, lead- 
ing to r = £a~^{cia + ci/3)[cosh(M;aA?/') — l]/2. However, 
the steric effect due to finite surfactant volume becomes 
important with increasing F. 

In this letter, we are interested in the equilibrium inter- 
face profile in the ID case. We require homogeneity of the 
surfactant chemical potential Ti>i = (SF/5ni)^ and the 
chemical potential difference Th = a^{dF/Sip)ni of the 
mixture. Some calculations give 



Pl = ln(ci/Zam) - 9ltp - ^iVi \n{(j)A(t>B) 

h = \n{(l)A/(l)B) - 2xV - C'VV - 5ici + K^,- 



(9) 
(10) 

The /lam stems from the amphiphilic free energy Fam = 
—T J drm InZam- Its ID form is 



Kin{z) 



dC 

2£ 



X{z + 



where X(z) = Ci{z) / Z am{z) . In the dilute limit Ci^ 
c\{z) is expressed in terms of ^'(-2) as 



ci{z) 



4iP{zf 



Ni/2 



Za.^{z)e 



(11) 



0, 



(12) 



In our model, if A''i is large and if ip^ is close to 1/2 far 
below the critical point, the first factor in the right hand 
side of eq. (12) can be large in the interface region where 
Ip ^Q. This leads to F ~ ci„^(l - A%Pl)-^^l'^ even for 
Zmn = 1- For our choice x = 3 in our figures, we have 
ipa = 0.429 and (1 - A.ipl)-'^ = 3.79. However, the first 
factor in eq. (12) tends to unity near the critical point. 

The grand potential f2 is given by the space-integral of 
w = /non — Ta~^{h'ip + vci). In equilibrium, O is mini- 
mized as a functional of ip and Ci under given boundary 
conditions. For an interface, u!{z) should tend to a com- 
mon limit Woo as 2: ^ ±00 (see the appendix) and the 
surface tension is expressed as 7 = / dz[uj{z) — lOoo]- Here 
it is convenient to introduce rj = a^fnon/T — y\Ci + Ci = 
a^oj/T + hip + ci. Use of eq. (9) yields 



V = {^+i')ln(pA + {^-ip)in(pB+NiCi 

-x^'-Jx^^? + lc(^')^ 



(13) 
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where ip' — dip/dz. The 7 can be expressed as 
T 



0.05 



7 



dz 



ri{z)-ria-h{^j{z)~ipa)-ci{z) + cic 



, (14) 



where h = (??a-'7,3-ciQ+ci;3)/(V'Q-'0/3) with 77^, rj/s, cia, 
and cif3 being the bulk values. The integrand in eq. (14) 
tends to zero far from the interface. We may then derive 
eq. (2) for small ci. Let tpoiz) — limci^oV'(^) be the 
interface profile without surfactant and Sip = ip—ipohe the 
deviation due to ci. Then r]{z)—r]a = Cd/dz[S'ipdipo/dz] + 
• • • up to first order in Stp and ci , because of the interface 
equation ln[{l + 2ipo) / {1 - 2ipo)] - 2x^o - C d'^ ipo / dz"^ = 
without surfactant. Hence a^Aj/T is nearly the integral 
of -ci{z) + cia - h{i/j — ipa) with h — [cip - Ciq)/A-0, 
leading to eq. (2). 

In all the figures to follow, we will set C = xa^, x — 
3, and I = 2.5a, where 70 = 0.497Ta-2 at ci = 0. In 
eq. (6) we have the critical density = 0.36ti;~^. The 
interface width ^ is about 5a 21. In fig. 1, we present 
ci(z) = a^ni(z), (70 —^)a?/T, and To? , where 51 = 4 and 
uia = 10. With increasing ciq,, the peak of ci(z) grows, 
r increases, and 7 decreases. The steric effect due to the 
surfactant volume fraction is apparent for A^i — 5, while it 
is not for A^i = 1. For A^i ^ 5, the first factor in eq. (12) 
is crucial for ciq, < 10^^ and the steric effect is relevant 
for larger ciq, so that eq. (2) holds only for very small 
ciQ;in fact, (70 — 7)/7T 1.3 at ciq — lO"**. See the 
argument below eq. (12). For A^i = 1, eq. (2) nicely holds 
for cia < 3 X 10""* and F oc cia at any ciq.. 

Ionic surfactant and counterions in a binary mix- 
ture. — Previously, we treated non-amphiphilic ions, in- 
cluding the solvation and image interactions [7,8]. Here 
the first ion species with density ni is a cationic surfac- 
tant. The second species with density n2 constitutes an- 
ionic counterions having no amphiphilic character. They 
are both monovalent with charges ±e. We assume that 
the counterions are very small and their volume fraction 
is negligible. Then the relations in eq. (3) still hold. 

For our complex ionic system, the total free energy den- 
sity / = /non + /ion IS the sum of /non in cq. (4) and the 
following new part. 



/ion = Tn2(lnc2 - .92 V") + 



(15) 



— a^n2. 



where n = rii + ri2. We define ci — a^ni and C2 
The dielectric constant is assumed to be of the linear 
form £{ip) = £c + Si'fpy where Sc and Ei are positive con- 
stants. The parameters gi in eq. (4) and 32 in eq. (15) 
can arise from the solvation (ion-dipole) interaction and 
can be very large in aqueous solutions [7,8]. The differ- 
ences A/i^^ = TgiAip are the so-called Gibbs transfer 
free energies per particle from phase a to phase f3 for 
ion species i [10]. For water-nitrobenzene mixtures at 
room temperatures, A/LtJ,^/T are of order 15 for mono- 
valent metallic ions. Recently, Sadakane et al. [11] found 
periodic structures in a near-critical mixture containing 
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Fig. 1: Results for nonionic surfactants for x = 3, gi = 4, and 
= 10: c\iyZ) for four values of c\a for A^'i = 5 (top) and 1 
(middle), while (70 — '))a? /T and Va^ versus cia for A^'i = 5 
and 1 (bottom). Here 70 — 7 — TP at small cia. 



strongly hydrophilic and hydrophobic ions, presumably, 
with gi ^ —g2 ~ 15. 

The electric potential <I> arises from the charge density, 

-V • e(V')V<I> = 47re(ni - 712). (16) 

The electric field E = —d^/dz around an interface reads 

47r^B 



-Eiz) = 

T ^ ' l + ei^(z)/ec j_oo 



dz'[m{z')-n2{z% (17) 



where the lower bound of the integration is pushed to — cx3. 
We define the Bjerrum length (.-q = /ScT at e = Ec- We 
may set = without loss of generality. 
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The image potential /iim acts on each ion [12, 13], which 
arises from inhomogeneous e. For moderate inhomogene- 
ity and in the ID case, it is expressed in the foUowing 
Cauchy integral [8], 



dz' e-2«l^-2'l d^{z') 



dz' 



(18) 



(19) 



where A represents the charge stength, 

and n = [Ame^ / ecT]^/'^ is the Debye wave number [12]. 
The factor g-^Kl^-z | ^j^g-j g^j-jggg from the screening 

of the image potential by the other ions, so the image 
interaction is weakened with increasing the ion density. 
We take n in k as the space-dependent ion density n{z) = 
ni{z) + n2{z). 

Again we assume homogeneity of the chemical poten- 
tials Ti/itot = SF/Sni, Tz^2tot = SF/dn2, and T/itot = 
a^SF/5ip. From F — J dr{f non + /ion) they are 

i^itot = 1^1+ e$/r + Aiim/r, 

1^2101 = lnC2 + 1 - 52-0 - e$/r + /iim/T, 
htot = - .g2C2 + ft-im/T, (20) 

where and h are given by eqs. (9) and (10). In taking 
the derivatives with respect to ni and n2, we neglect the 
n-dependence of in eq. (18). The hi^{z) in eq. (20) 
arises from the image interaction and is obtained by the 
right hand side of eq. (18) if tp is replaced by rt = ni -I- 712. 

Since the counterion volume fraction is neglected, C2 
around an interface is written as 



C2a 



exp 



e 1 



(21) 



In our previous work [8], we examined the image factor 
-Fima(z) = exp[— /iim(2)/T] in various cases. Let us con- 
sider it in the a region with Zin ~ z > ^, where Zin is the 
interface position and ^ is the interface thickness. Then, 

Fi„,a(2) = exp[-i?e-2-=(-.n-.)/(^.^^ _ ^)]^ (22) 

where D = iBSiAt/j/Asc- Here D ^ £b ^ Aa for ei ~ Sc 
and Alp 1. The image interaction can be crucial near an 
interface for ^ < Z? < under which the image factor 
serves to repel the ions in the region ^ < Zjn — z < (2kq)~^ 
in the a region. This ion depletion was used to explain 
an increase of 7 of water-air interfaces with salt [12, 13]. 
However, its elect on the ionic surfactants is diminished 
when the amphiphilic interaction is strong or for WaAip ^ 
1. These aspects will be illustrated in figs. 5 and 6 for 
^ = 4 and 10. 

From the charge neutrality in the bulk regions a and /3, 
we require cia — C2a and cip = C2/3. Then we find 

ln(ci„/ci^) = (.gi + g2)A^/2 + St, (23) 

e(*a-$/3)/T=(gi-52)AV'/2 + 5t, (24) 
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Fig. 2: Profiles for mixtures with cationic surfactant and an- 
ionic counterions with A^'i = 5, Wa = 12, and A = 4. Top: 
A'^ici (bold line), <f)A, and (j>B- Middle: ci and C2. Bottom: 



e^/T. Here gi = 4, g2 = 10, and ciq 



10"'' (left), while 



gi — — 32 = 8 and cia = 3.6 x 10~* (right). The counte- 
rion distribution has a peak in the phase a (left) or [3 (right) 
depending on g2. 



where — $^ is called the Galvani potential difference. 
The St stems from the surfactant volume fraction NiCi, 



St 



In 



(1 - 7ViCi<,)2 - AlPl 



4^1 



xAci, (25) 



with Aci = Cia — cip. Here 5't ^ as iVi — > 0. Note that 
$(z) changes on the scale of the Debye screening length, 
in the phase a and in the phase f3. As a result, $ 



changes from to $^3 on the spatial scale of 

We consider the grand potential density u — f — 
Ta''^{htotip + i^itotCi+V2totC2)- With the aid of eq. (20) we 
have u) = Ta^^{r] - /itotV" - c) + £| V$|^/87r - p$, where 77 
is given by eq. (13) and c — ci +C2. From lu{oo) = uj{—oo) 
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Fig. 3: 7a^/r, ^\a? jT , and Fa^ as functions of ci^ with A'^i — 
5, X — 3, and uia = 12, where 7 is calculated from eq. (26) and 
71 is the first term on its right hand side. The curves change 
on a scale of 10"'' for gi = 4 and gi = 10 (top) and on a scale 
of 10~* for (ji = —g2 = 8 (bottom). 



(see the appendix), we obtain 



7 = ^ / 



77(2) - htoti'iz) - c(z) - Aa 



y dz£(^)|V$| 



(26) 



where Aa = ?/q - /itotV^a - Cq and /itot = - »y/3 - Cq + 
0^)1 IS.1I) with Cq and being the bulk values of c. From 
this expression, eq. (2) follows at small n = ni+n2 = a^^c 
if we use the argument below eq. (14). 

We present some numerical results with A^i = 5 and 
El = O.SSc- The dielectric constant of the phase a is twice 
larger than that of the phase /3 at x = 3. We set ^ = 4 
except in the left panel of fig. 6 where A — 10. 
(i) In fig. 2, we show the volume fractions (j)A, <j)B, and 
NiCi — \ ~ 4>A — 4>B (top), ci and C2 (middle), and e$/T 
with $Q, — (bottom). In the left plates, the counte- 
rions are more hydrophilic than the cationic surfactant, 
where 51 = 4 and 52 = 10 leading to F = 0.124a~^ and 
7 = 0.317Ta"2 at ciq = 10"^. In the right plates, the 
surfactant cations are hydrophilic and the counterions are 
hydrophobic, where gi = —52 = 8 leading to F = 0.155a~^ 
and 7 = 0.159Ta-2 at ci^ = 3.6 x lO"". The distribution 
of the ionic surfactant ci is narrower than that of the coun- 
terions C2. This gives rise to a peak of $, at which the right 
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Fig. 4: ja^ /T, ■jio^/T, and Fa^ as functions of g2 with x — 3, 
Wa = 12, and gi = 8, where 71 is the first term on the right 
hand side of eq. (26). The counterions are hydrophobic for 
g2 < and hydrophilic for (72 > 0. 



hand side of eq. (17) vanishes. 

(ii) In fig. 3, we show 7, 71, and F as functions of Cia at 
Wa — 12, where 71 is the first term on the right hand side 
of cq. (26). The Lifshits point cia = a^riL = 2.5 x 10^^ 
for C = 9a^ is marked by an arrow (top). For small cia 
eq. (2) holds, but the last electrostatic term in eq. (26) is 
not negligible. For a pair of hydrophilic and hydrophobic 
ions, a large electric double layer is formed at an inter- 
face [8], leading to a large F even at small ciq. 

(iii) In fig. 4, we show 7, 71, and F versus g2 at Wg. — 12 
and (71 = 8 to demonstrate the above trend, where F and 
70 — 7 are markedly enhanced for negative g2 ■ 

(iv) In figs. 5 and 6, we set Wg, ~ 15, where F is ap- 
preciable even at cia ~ 10^^. In fig. 5, where gi — 
32 = 10, we have F = O.OOgOa^^ and 7 = 0.487T/a2, 
while 71 = 0.491T/a2 70 - TV. There are virtually 
no ions in the region /3. The distribution of C2 in the 
region a is broad changing on the scale of where 
Ka ~ 0.046/a. Note that ^{z) is constant in the region 
f3 displayed in fig. 5. However, it should tend to ^i3{= 
here) for z — Zin S> lO^a. 

(v) In the left panel of fig. 6, we increase A in eq.(19) to 
10 with the other parameters being the same as in fig. 5. 
This is because the effect of the image interaction is rather 
weak at A = 4 [8], where D in eq. (22) is of order a. See 
the discussion around eq. (22). With A — 10, the image 
interaction is amplified and F is decreased to 0.0049a~^ 
(55% of the value in fig. 5) with 7 = 0.491T/a2. See 
fig. 10 of our previous paper [8], where we set A = 10 to 
realize strong ion depletion in the more polar phase. How- 
ever, the normalized potential e^/T is also amplified with 
increasing A from eq. (17) and the effect is very compli- 
cated. 

(vi) In the right panel of fig. 6, we consider the case of a 
hydrophilic and hydrophobic ion pair with gi = —52 = 8. 
The other parameters are the same as in fig. 5. Then F 
is incresaed to 0.071a~^ (8 times larger than the value in 
fig. 5) with 7 = 0.403T/a2. 
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Fig. 5: Ci and C2 (left) and e$/r (right) with = 15, yl = 4, 
and ffi = <72 = 10, where c\a = 10~^ is very small due to large 
Wa- The counterion distribution is wider than that of the ionic 
surfactant. There is no appreciable variation of $ in the /3 
region shown here, where the ion densities are nearly zero. 
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Fig. 6: ci and C2 for — 15 and ci„ — 10~^. Left: A = 10, 
where the electrostatic and image interactions are amplified 
from eqs. (17) and (18). Right: gi = —g2 = 8 leading to 
enhanced adsorption. The other parameters in each panel are 
the same as in fig. 5. 



Summary. — We have presented the continuum mod- 
els of surfactants including the amphiphilic interaction ex- 
plained around eqs. (4)- (8). They reasonably describe the 
adsorption of nonionic and ionic surfactants, though our 
numerical examples are still fragmentary. For ionic sur- 
factants, included also are the electrostatic, solvation, and 
image interactions. The adsorption is extremely sensitive 
to Ca/T = Wi^Atlj and Ni. It is decreased as ^ ~ ^b/o is 
increased for a pair of hydrophilic ions from the left panels 
of figs. 5 and 6. It is also enhanced for hydrophilic and 
hydrophobic ion pairs, as shown in figs. 3 and 4 and the 
right panel of fig. 6. 

Though we have set x = 3 in all the figures here, the 
changeover and the phase transition with decreasing x at 
large Wa. and gi should be intriguing [8]. The phase tran- 
sition of binary mixtures with salt can be very complex at 
large gi [7,8,11]. We should also study dynamics including 
the amphiphilic interaction in higher space dimensions. 

Appendix. — We consider a ID eqilibrium interface 
profile, ■!/; — iIj{z), ci — Ci(z), and C2 = C2{z), for our ionic 
surfactant system. Use of eqs. (9), (10), and (18) gives 



dz 



T 



T T' 



where c = ci + C2, = dip/dz, = dcx/dz, and 



T 5il) ' 



^Kint — 



T SCK ' 



(A.l) 



(A.2) 



with K = 1,2. Here Fint = / dr[—Tni\ia Z^m + "-A'im] 
is the sum of the amphiphilic and image free energies, 

so hint = /lam + him and Z^iint + In Zam = i^2int — Mini 

(see the discussion below eq. (20)). The Fint is invariant 
with respect to a shift of the interface position or with 
respect to z ^ z — in ip and ck- Then J dzlhiatip' + 
i^iintc'i + i^2intC2] = and integration of eq. (A.l) gives 
w(cxd) = a;( — oo). 
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